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ABSTRACT 



The purpose of this thesis is to research the 
availability of targets moving on or near the earth's 
surface when viewed by an orbiting satellite. A discussion 
of basic orbital mechanics is presented as well as a 
development of a suitable coordinate system. An analysis of 
non-linear observability is then provided. Lastly an 
observer is designed and sucessfully simulated. 
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I . INTRODUCTION 



This thesis deals with satellite tracking of relatively 
low altitude targets. The long term goal of this project is 
to determine an effective general -coverage satellite orbital 
pattern. The short term goal of this thesis is to explore 
the possible use of various coordinate systems and to design 
a satellite observer. A target such as a slow moving 
aircraft is assumed for the basic development. 

Chapter 2 contains an overview of basic satellite 
mechanics and dynamics as well as a brief discussion of 
satellite detection and tracking equipment. Equations 
regarding orbital motion will also be introduced. Also 
chapter 2 provides an overview of the entire thesis. 

The theory and mathematics of satellite and target 
motion for a specific case are developed in Chapter 3. The 
geometry of motion on a spherical surface is detailed in 
that chapter. Specifics such as geometric swath width, 
line-of -sight and precession are furnished. This leads into 
an in-deapth examination of coordinate systems. Various 
attempts to find a suitable set of reference axes for 
analytic studies of the tracking problem are presented. It 
is desired to achieve a target and observation model that 
has a convenient basis for analysis such as either linear or 



bilinear in state space- Several trial coordinate systems 
were developed mathematically but found unsuitable. One 
trial system is considered satisfactory under certain 
assumptions and is, therefore, completely developed. 

Chapter 4 presents a discourse on observability in 
general, and non-linear observability in particular. The 
concept of connectedness and univalence is discussed. The 
chosen system is then analyzed for observability. 

The topic of observers is presented in Chapter 5. A 
background on basic obsevers is provided. Finally an 
observer is designed and simulated by computer for the 
prefered system. 

Chapter 6 offers conclusions and recommendations for 
further studies. 
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II- THE SATELLITE IN ORBIT 



Before an understanding of satellite tracking can be 
achieved it is essential to have an understanding of basic 
satellite dynamics and satellite sensors. The concepts and 
equations that govern the motion of a satellite in orbit 
around the earth are based on the physical laws put forth by 
kepler and Newton- These laws of motion apply to artificial 
satellites as well as planets and moons. The physical 
geometries and forces are the same. The following 
explanations and derivations form a basis for orbital motion 
of a satellite. 

A. THE SHAPE OF THE EARTH 

In reality the earth is neither exactly round nor 
exactly a rigid body. The constituents can be molded and 
deformed to a certain degree. The earth's rotation causes 
it's own materials to feel a force that pulls tangentially 
into space. Fortunately the earth's materials also undergo 
a constant centripetal (center seeking) acceleration to keep 
them in place. This centripetal acceleration, which is 
directed towards the center of curvature, is provided by 
earth's gravity. 

The gravitational force produced by the earth is 
directed towards the center of the earth. However, the 
momentum produced by the earth's rotation pulls at the 



earth's materials. Each particle of earth has both these 
forces acting on it. The earth's surface at the equator is 
actually moving faster around the rotation axis than the 
surface material at or nearer the poles. Figure 2.1 
illustrates this. Therefore the materials on the surface at 
the equator feel this tangential pull the most. Analysis 
shows that the result is an earth shaped like an oblate 
el 1 ipsoid . 

An oblate ellipsoid can be described as a sphere which 
has been compressed along the polar axis and therefore 
bulges at the equator. A cross sectional slice through the 
poles yields an ellipse. A slice along the equatorial plane 
yields a circle. 

This oblateness (measure of the earth's flattening) of 
the earth is only slight. The equatorial diameter is 
calculated to be 12,757 km whereas the polar diameter is 
calculated to be 12,714 km. This is only a 43 km difference 
which is about 1 part in 297. 

B. GREAT AND SMALL CIRCLES 

The earth's equatorial bulge can be temporarily 
disregarded in order to consider the earth as a sphere 
turning beneath an orbiting satellite. The intersection of 
a plane passed through the center of the sphere and the 
sphere is the largest circle that can be drawn on the 
sphere's surface. This is known as a great circle. The 
shortest distance between any two points on the surface of 
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the earth is an arc of a great circle- The earth's equator 
is an example of a great circle- There are an infinite 
number of approximate great circles on the earth's surface 
since there are an infinite number of planes that can be 
passed through the earth's center. 

There is a series of 24 special great circles that pass 
through the earth's poles and are evenly spaced from each 
other. These are called meridians and they intersect with 
the equator at right angles- Meridians join at both poles. 
Forming other right angles to the Meridians are parallels. 
Parallels are small circles created by passing planes 
through the earth parallel to the equatorial great circle. 
See Figure 2.2. 

In order to determine precise locations on the earth's 
surface, latitudes and longitudes are utilized. The 
longitude of a specific place refers to the arc (measured in 
degrees) of a parallel between that place and the prime 
meridian (which passes through Greenwich, England for 
reasons of history). Longitudes run east and west. 
Latitudes, however, run north and south. The latitude of a 
specific place may be defined as the arc (in degrees) of a 
meridian between that place and the equator. Figure 2.3 
shows an example. When the earth is considered as having an 
ellipse as a cross section (which it has) instead of a 
circular one, the length of a degree of latitude is slightly 
greater at the poles than at the equator. 
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Figure 2-2 Meridians anci Parallels 
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C. ORBITAL MOTION 

A presentation concerning orbital motion is basic to 
satellite tracking. Gravitational attraction and momentum 
combine to keep a satellite in perpetual orbit above the 
earth (in the absence of friction and external celestial 
attractions) . 

1 . Gravity 

The understanding of gravity and orbital motion can 
be traced back to Johannes Kepler (1571-1630) and Isaac 
Newton (1642-1727). Kepler is known predominantly for his 
three laws of planetary motion. The laws are as follows; 

(1) Each planet moves about the sun in an orbit that is an 
ellipse, with the sun at one of the foci of the ellipse. 

(2) The straight line joining a planet and the sun sweeps 
out equal areas in space in equal intervals of time. 

(3) The squares of the sidereal periods of the planets are 
in direct proportion to the cubes of the semimajor axes of 
their orbits. 

The above three laws apply equally well for 
artificial satellites as they do for planets. Newton 
restated and clarified Kepler's laws. Newton was an 
advocate of rigorous proofs whereas Kepler preferred to 
state empirical laws based on observations. 

Kepler's first law stated mathematically is 

f = mv 2 /r (2.1) 
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where: f =centr ipetal force needed for circular orbit, 

m=mass of planet (or satellite) in orbit, 
v = velocity of orbiting body, 

r = distance between sun and planet (or satellite and 
earth) , 

Kepler's third law restated mathematically is 

(m f + m^p a = 4ir a a /G (2.2) 

where ml and m2 refer to the masses of two bodies that 
revolve mutually about each other, 
p = period of revolution, 
a = semimajor axis of relative orbit, 

. i 

G = universal gravity constant ( 6.67 x 10 ) 

Newton also explained the reasons behind kepler's 
original observations. That is, Newton supplied the laws of 
motion which are at the root of classical mechanics. 

These laws are as follows; 

(1) Every object remains at rest or in uniform motion 
unless an external force acts upon it. 

(2) The product of the mass of an object and its 
acceleration vary directly as the resultant force, and 
the change in motion takes place in the direction of 
that force. 

(3) For every action there is an equal (in magnitude) and 
opposite (in direction) reaction. 

In the case of orbiting planets or satellites, 

gravity is the unseen force referred to in Newton's first 

law. Newton also postulated that every particle of matter 

in the universe attracts every other particle. The force of 
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this attraction is proportional to the product of their 
masses and inversely proportional to the square of the 
distance which separates them. 

In equation form this is 

F = Gm, m z /d J (2.3) 

where d is the distance between the center of mass 1 and 
mass 2. 

Orbital motion can be understood in terms of 
Newton's laws previously stated. The critical factors in 
putting a satellite in orbit (by the traditional launch 
method) are speed and direction of movement at burnout. 
Burnout is when the rocket engine shuts off and the 
satellite behaves as an astronomical object. 

As a satellite follows it's orbital path it 
continuously falls toward the earth due to the earth's 
gravitational pull. However, the satellite's momentum 
prevents it from really being pulled any nearer to the 
earth. The satellite's orbit is the result of two main 
forces. The momentum of the satellite is a measure of its 
state of motion. The inertia of the satellite (recall 
Newtons first law) is that property that causes the 
satellite to resist acceleration and travel in a straight 
line. For the satellite to move in a circular path rather 
than in a straight line, it must continually suffer an 
acceleration toward the center of the circle. This 
acceleration is centripetal acceleration. The central force 
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that produces the centripetal acceleration is that 
gravitational attraction between the satellite and the 
earth. Figure 2.4 illustrates this in vector form. 

The period of an artificial satellite is directly 
related to the size of it's orbit; i.e, the closer the 
satellite is to the earth the faster it travels. 

2 . The Center of Mass 

The center of mass can be defined as that point 
within a system that either remains fixed or moves as if the 
entire mass of the system were concentrated at that point. 
For a number of particles the center of mass is defined as 

r = (l/M T<rb . L ).Zm- r. (2.4) 

An example of this for a two body system is 

x c ,rv = (2.5) 

m“ + ’ m* 

Mass 1 is a distance x(l) from an arbitrary origin, 
and mass 2 is a distance x(2) from the origin. X(cm) 
represents the location of the center of mass. A solid 
object can be thought of as a collection of a great many 
particles. In this case equation 2.4 applies. 

Very often it is useful in certain calculations to 
treat a particular body as if its entire mass is 
concentrated at one point. That one point is the center of 
mass • 
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Figure 5.4 Centripetal deceleration 
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It can be shown [Re£. 1] that a large sphere 



attracts other bodies as though the sphere's mass is all 
located at the center- This holds true as long as the 
sphere is of uniform density or is made up of concentric 
shells each being of uniform density. The earth is often 
approximated as being a sphere made up of uniformly dense 
shells [Ref. 23 . 

3 . Equations of Motion 

A satellite in orbit around the earth will feel by 
far the strongest gravitational pull from the earth itself. 
Other astronomical bodies, such as the m’oon or sun, do exert 
force disturbances but are too far away to exert a very 
strong pull on the satellite. Assuming the gravitational 
acceleration of the satellite is due only to the earth's 
gravitational attraction, g can be defined as 

g = Gm^/r 2 (2.6) 

The value of g does vary with respect to altitude 
and latitude (since the earth is not really an exact 
sphere). Tables 1 and 2 depict the variations of g. The 
values of g are close enough to allow an approximation of 
the earth's shape as spherical for most applications. 

An orbiting satellite is considered to be in the 
earth's gravitational field. The vector symbol g refers to 
this field and is defined as 
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TABLE 1 



VARIATION OF g WITH LATITUDE AT SEALEVEL 





Latitude 


g (m/s i ) 




0” (equator) 


i 

9.78039 


1 


10* 


9.78195 




20° 


9.78641 




30* 


9.79329 




40° 


9.80171 




50* 


9.81071 




60* 


9.81918 




70* 


9.82608 




80* 


9.83059 


( 

i 


90* 


9.83217 



‘ 


TABLE 2 


- 


VARIATION OF g WITH ALTITUDE 


AT 45° LATITUDE 


r 

Altitude < km ) 




i 

g(m/s 2 ) 

| 


0 




9.806 


1 




9.803 


4 




9.794 


8 




9.782 


16 




9.757 


32 




9.710 


100 




9.600 


500 




8.530 


1000 




7.410 


380000 




0.003 

J 
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g 



F/m 



< 2 . 7 ) 



When on object is a distance h above the earth's 
surface it has potential energy defined as 

U = - ( -mg) h = mgh (2.8) 

Gravity is a conservative force pointing to the 
earth's center and has the value (-mg). The gravitational 
force exerted on an object in this field can be derived from 
the potential energy equation as 

F = -du = d f-GMml = -GMm (2.9) 

dr dr L r "j r 2 

A more general form for potential energy invoking 
the universal constant G is 

U = -GMm/r (2.10) 

Potential energy can be converted into kenetic 
energy [Ref. 3]. As an object falls to earth it loses its 
potential energy as its height decreases but gains kenetic 
energy as its velocity increases. Kenetic energy is 
mathematically defined as 

K = GMm/2r = mv 5 /2 (2.11) 

The total mechanical energy is 

E = K «■ U (2.12) 

Orbital motion can be considered as motion in a 
plane. For one full revolution around the earth the 
satellite remains in a plane as it traces out an ellipse or 
circle . 
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The velocity, period and altitude of a satellite are 
all inter-related. The period of a satellite [Ref. 23 can 
be defined as 

7T= 2n (R + h)/v , <2.13) 

where v is the velocity of the satellite. 

The centripetal force, is often defined as 

F = m 3 v i /(R + h) <2.14) 

where m is the mass of the satellite. 

It is obvious to the casual observer that the 
centripetal force is also in part a function of the 
satellite's altitude and velocity. Polar coordinates are 
used in figure 2.5 to illustrate the elliptical orbit 
variables . 

4 • Uniform Circular Motion 

Sometimes a satellite's elliptical orbit is that 
perfect ellipse the circle. If the satellite's orbital 
velocity remains unchanged in a circular path then the 
satellite is moving with uniform circular motion. 

Since the orbit of a two body system <earth and 
satellite) is often in a plane it is possible to use polar 
coordinates to develop basic equations. Figure 2.6 shows 
the relation between polar and rectangular coordinates. 
Velocity in polar coordinates for circular motion can be 
expressed as 
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Figure 2.5 Elliptical Orbit 




Figure 2.6 Polar and Rectangular Coordinates 
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V 



(2. 15) 



= r 0 0 

where 0 refers to motion in the theta (0) direction. 
Acceleration is found by taking the derivative of the 
velocity 

a = -r0 2 r + r08 C2.16) 

Energy can be expressed as 

E = K + U = m(r J 0 J )/2 - GMm/r C2.17) 

and angular momentum as 

L = mr J 0 C2.18) 

Uniform circular motion is easier to deal with than 
non-uniform circular motion. There is only an angular 
component of velocity present in the equation. It is often 
useful to simplify the satellite dynamics in order to more 
easily manipulate equations. 

5 . Equations of Motion in Polar Coordinates 

Consider the special case of a satellite in a 
circular orbit. This can be treated as motion in a plane. 
Polar coordinates are especially usefull for this case. The 
two variables that define motion and position are the radial 
coordinate, r, and the angular coordinate, 0, conversion 
between rectangular coordinates and polar coordinates is as 
follows . 

r = \) X 1 + Y 1 X = rcos C 8 ) 

0 = arctanCY/X) Y = rsin(0) 

Figure 2.7 illustrates this. 
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Polar coordinates are often used any time 
curvilinear motion occurs- In Figure 2.8 basic curvilinear 
motion is shown. The differential length dr has components 
in both the r and 9 directions. The vector r can be 
expressed as rn where n is a unit vector in the radial 
direction . 

Velocity can be defined as; 

v - d(r) = d(rn r ) = d(r^ n r ) + rdn r (2.19) 

dt ~dt dt dt 



Rectangular coordinates are used once again to achieve an 
expression for the time derivative of the unit vector n and 

He ; . 

n_ r = cos<9>i ♦ sin(9>} = -sin<9)i + cos(9)] 

n r = -sin<9)i9 ♦ cos(9)j9 = <-sin<9)i ♦ cos(9>j>9 = 8n e 
ho = -cos(9)i9 - sin(9)j9 = -(cos(9)i + sin<9)j)9 = -9n r 

Velocity can now be expressed as; 

V = rn ♦ r9n„ 

-r -9 

Following the same rules of differentiation acceleration can 
be derived. 



a = V = rn * r9n + r9n, +r9n,. -r9'n 

— - -r -3 -e -o 

Simplifying , 

a = <r - r9 1 >n ♦ (2r9 * r9)n_ (2.20) 

— - r 

For circular motion the radial component, r, remains 
constant. Only 9 varies. Therefore; 

( 2 . 21 ) 



V = r9 ^ e 



a = -r9 2 n ♦ r9n 



— 6 



< 2 . 22 ) 
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Fig ure 2. 3 Curv ill near 
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The acceleration term directed in the negative radial 
direction is the center seeking acceleration. 

Figure 2.9 depicts the directions for the force vectors. 
The sum of forces is calculated as follows; 

ZF = ma = m(r-r9 * ) n C2.23) 

- r* — r — r 

£F . = ma = m(r9 + 2r9)n C2.24) 

-e -e - e 

Using differential equations and state form motion can be 
represented as in equations 2.27 and 2.28. 







dr/dt = V r 


(2.25) 






de/dt 


(2.26) 


dV 


/dt 


= V* /r - k/r «■ U - 6,- 


(2.27) 


d- 


/dt 


= -2V r i ^/r +OJ -6 ft )/r 


(2.28) 



where U f = F r /m and = F^/m. Both variables represent 
thrusts. The symbols 6 r and 6^ are used to represent 
disturbances which are suumed to be zero in this case. 

D. SATELLITE SENSORS 

Various sensors on board the satellite are employed to 
keep the satellite itself on track, and to detect and track 
a possible target. The satellite guidance system uses 
sensors for measuring certain vehicle dynamic variables. 
Such variables include satellite acceleration, velocity, 
position and angular velocity. Guidance systems generally 
determine these variables, compare the received information 
with desired parameters, and generate correction commands, 
various types of sensors include inertial y optical and 
radio sensors. The satellite may make use of gyroscopes to 
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establish a set of reference axes on board. Active thrust 
control is used to align the body axes with the reference 
axes. Once the satellite has established the craft's 
position and attitude active control components are 
activated. These components respond either to a telecommand 
from a ground station or to the satellite's on-board 
computer in order to correct the satellites's status. In 
some applications where perturbations are likely, feedback 
control is added to the active control of a satellite. In 
some cases optical sensors are used in conjunction with 
inertial equipment. Optical sensors include such devices as 
sun finders, sun sensors, planet sensors, celestial 
trackers, and horizon scanners. The gyroscopes and 
accelerometers are examples of inertial sensors. In certain 
applications the inertial sensors provide short-term 
stability and optical sensors provide long-term stability. 

The gyroscope is an inertial guidance system that is 
often used for ships, aircraft and spacecraft. A gyroscope 
is a device that possesses a high rate of spin about an axis 
of symmetry that has freedom of angular rotation. Examples 
of gyroscopes are the rate gyro and the integrating gyro. 

Surveilance satellites normally have an entirely 
different set of sensors for target acquisition. Very often 
this type of satellite carries several types of information 
gathering devices. High resolution optical devices, radar 
scanners and infrared capacity can all be included in a 
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satellite payload. Passive systems are sensitive receivers 
which normally provide data on bearing and bearing rate. 
These systems have the advantage that they do not alert the 
target of their actions. Optical techniques normally 
involve lenses and cameras. This can be very precise, 
especially at low altitudes. 

Radar is an example of an active sensor. It bounces 
radio waves off a chosen target. This is excellent for 
determining both bearing and range. Both optical and radar 

methods can be degraded by poor atmospheric conditions. For 

% 

the radar tracking technique noise is always the main 
limitation. 

Inverse scattering methods are being developed to enable 
polarimetric radar to obtain better target information. 

This procedure involves illuminating the target with 
polarized waves and observing the amplitudes and phases of a 
set scattered waves. 

Another type of sensor is the synthetic aperture radar 
system. This is based on holography methods. Holography 
has the ability to record three dimensional pictures and 
focus sharply on both the near field and the far field at 
the same time. This involves recording a wave interference 
pattern. A microwave generator is used to provide a 
constant frequency microwave signal as well as a reference 
wave. This is considered to be a highly precise and 
accurate sensor . 
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Passive sensors are normally infrared sensors- There 
are two major techniques used in enfrared detectors- The 
two main types are called photon detectors and thermal 
detectors • 

In photon detectors, the technique used involves a 
photon of infrared radiation being absorbed by a semi- 
conductor electron which raises its energy level into a 
conduction band- Photons with less energy than the band gap 
produce bo signal and are thereby effectively filtered out. 

Thermal detectors function be sensing the temperature 
change resulting from absorption of infrared radiation by a 
suitable element- Usually, this absorbing element has some 
temperature sensitive electrical property such as 
resistivity so that the temperature change is sensed 
electrically - 

One device often used in passive detectors is the 
radiometer- This device is a broadband, dual frequency , 
low noise, solid state, remote controlled mechansim. 
Radiometers are used to obtain high resolution imagery from 
low altitude satellite's. 

A promising type of detector is the Silicide Schottky 
diode based infrared camera- These are easy to manufacture 
and have excellent performance- The camera focal plane 
consists of a two dimensional array of metal electrodes 
fabricated on a silicon substrate- The focal plane is back 
illuminated- When an infrared picture is observed, hot 
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carriers are emitted from the focal plane electrodes. An 
electronix image of the scene is formed by accumulation of 
these carriers on a pixel by pixel basis. 

Another' technique involves the Charge Injection Devices 
(CID) . These devices are surface charge devices that 
collect photon generated charges and store them in MOS 
capacitors . 

Other types of infrared detectors include the Lead 
Sulfide (PhS) Detectors, Lead Selenide Detectors (PbSe) , 
Thermister Infrared Detectors and Indium Antimode Charge 
Injection Devices . 

One tracking technique available to surveillance 
satellites involves the use of the doppler effect. The 
satellite emits a signal and then receives the signal's 
return after it reflects off a given target. The change in 
frequency as a result of the relative motion between 
satellite and target leads to a calculation of the targets 
velocity. Ref 3 defines the return frequency, ^ , detected 
by the satellite as 



V( l-u/c)/ l-(u/c> 



(2.29) 



where 

> = frequency detected if both were at rest, 

u = relative separation speeds, 
c = speed of light. 
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Artificial satellites are often subject to a variety of 
disturbing forces. These include atmospheric drag, 
variation of atmospheric density, solar radiation pressure, 
surface charge drag, meteorite impacts, lunar or solar 
gravity caused perturbations and possible encounters with 
hostile killer satellites. These potential disturbances can 
act seriously to impede the position as well as the attitude 
stability of a satellite. Fortunately, there are passive 
and active stabilization systems to enable recovery from 
most perturbations . 
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III. SATELLITE COVERAGE AND COORDINATE SYSTEM S 



This chapter includes satellite coverage, spherical 
geometry and coordinate systems. A coordinate system which 
enables analytical calculations of satellite observer and 
tracker is desired for at least a single case. Several 
systems are explored in an attempt to find a suitable set of 
reference axes to establish a base for more complicated 
analysis by computer. A simple case is introduced and an 
appropriate coordinate system is presented. 

A. G’ENERAL SATELLITE COVERAGE 

Basic information on satellite coverage is discussed 
here. An understanding of geometric swath width and 
precession is required for a detailed study of satellite 
tracking . 

1 • Geo metr i c Swath Width 

As a satellite travels around the earth it covers a 
certain amount of surface area. This area is in the form of 
a spherical cap (figure 3.1). As the satellite's altitude 
above the earth increases the area of the spherical cap 
increases. The area of the spherical cap [Ref. 5] can be 
calculated as 

A = 2ttR 2 ( l-sin( 90 ° - 0 ) ) (3.1) 
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where R is the radius of the earth and 9 is the angle shown 
in Figure 3.2. 

It appears from this that a high altitude orbit is 
beneficial. However, it must be considered that as 
satellite altitude increases, satellite sensor accuracy 
decreases . 

If a target is at point 'a' (as in Figure 3.2) it 
can be viewed by the satellite as long as the satellite is 
between points 1 and 2 on its orbit. This portion of the 
orbit keeps the satellite above the horizon with respect to 
the target. The amount of time the satellite is above the 
horizon is calculated [Ref. 53 to be 

t = 29 /w (3.2) 

where w is the angular velocity of the satellite. 

As the satellite progresses in its orbit the spherical 
cap of coverage moves with it. As it moves, the cap traces 
out a ribbon around the earth. The width of the ribbon 
(Figure 3.3) is the geometric swath width (GSW) . 

Figure 3.4 illustrates the GSW as the very dark arc on 
the earth's surface. The dotted line from point 1 to 
point 2 is tangent to the earth at the surface location 
directly beneath the satellite. R is the radius of the 
earth and h is the satellite altitude. The arc, a(l), in 
Figure 3.5 is one half the dark arc in Figure 3.4. Simple 
geometry proves that arc a<l) = R9 . Therefore the 
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i-igure 1 



Spherical Cao 




F i gure 



View of the Target 





Fig ure 



4 



Geometric Swath Width 



GSW = 29R. Figure 3.6 depicts a right triangle. From the 
rules of trigonometry cosine(Q) = R/(R+h). Therefore, 

9 = arccosine < R/ ( R+h ) ) 

It follows by substitution that 

GSW = 2R ( arccosine (R/(R+h)>) (3.3) 

2 . G round Track and Coverage 

As a satellite travels in its orbit the distance it 
moves in one full trip around the earth is 

D = 2tt ( R + h ) (3.4) 

where R is the radius of the earth and h is the satellite's 
height above the surface. The orbital speed is calculated 
[Ref. 5] to be 

4.2685<R/(R+h) ) km/sec (3.5) 

The satellite also has a ground speed and a ground 
track. The satellite's position can be projected on to the 
earth's surface. This point on the earth's surface directly 
below the satellite is called the sub-satellite point (SSP). 
As the satellite moves its SSP forms a ground track which 
can appear as simple as a circle about the earth or a mere 
point on the surface or they can appear very complicated. 

The ground speed due to the satellite's motion alone is 

S m s 7.9053(R/ (R + h) )km/sec (3.6) 

However, the earth does rotate at about 0.267472 
km/sec at the equator. The actual ground speed depends on 
the inclination of the orbit plane to the equatorial plane. 
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Figure 3.6 Right Triangle 



For all cases the ground speed is less than the actual 



orbital speed. As the ground track is formed so is the 
ribbon whose width is the GSW. Because of the earth 
rotation beneath the satellite the swaths may overlap at 
some places yet never cover other places. If the earth's 
rotation is ignored and its shape is assumed spherical a 
satellite staying in one plane traces the same circular 
swath on every orbit. 





If the earth's equatorial bulge is taken into 


account 


the plane of the orbit precesses about the earth's 


axis of 


rotation. This precession is caused by the non 


central 


gravity force field in which the satellite travels. 



It is measured by the rate of drift of longitude on the 



surface 


of the earth. 

For an elliptical orbit the rate of precession is 

ft = 9 . 95 ( R/a ) ^ ( cosine (i)/(l-e a ) 2 (3.7) 


where ft 


= orbital plane precession rate in degrees/day 


R 


= mean radius of the earth 


a 


= semimajor axis 


e 


= eccentricity of ellipse 


i 


= angle of inclination of orbit 

For the special case of a circular orbit the 



precession rate is 

ft = 9 • 95 < R/ ( R+h ) ) * cosine(i) C3.8) 

If the effect of this precession is ignored, the 
motion of a satellite lies in a plane passing through the 
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center of the earth. A satellite making one revolution per 
day at 0° inclination to the equator is stationary over a 
point on the equator. In this case a stationary spherical 
cap on the earth's surface is the only area covered. 
Commercial communications satellites are usually 
geostationary. This is a special case of the synchronous 
orbit <24 hour orbit) . If a synchronous orbit is inclined 
to the equator the ground track is a figure eight. 

3 . The Spherical Triangle 

At this time all the simplifying assumptions are 
invoked to establish a base for further analysis. The 
satellite is traveling in a circular orbit above the 
equator. Obviously, transformations to other orbits are 
readily available. The SSP is always on the equator. There 
is a target traveling on a steady course at a constant 
velocity and therefore following a great circle heading. 

This is not an unreasonable assumption. It is normal for a 
ship or an aircraft to keep to a great circle heading at a 
most efficient cruising speed and altitude. The target is 
assumed to be on the earth's surface or very near to the 
surface. Figure 3.7 illustrates the paths for both the 
surveilance satellite and the target. The target track is 
inclined i* to the satellite track. It is assumed that both 
can be considered to be moving on the surface of the same 
sphere with transformations readily available to the actual 
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F i gure 3. 7 



Satellite and Target Track 



satellite position. The earth is assumed to be a perfect 
sphere and is approximated by a point mass at the center. 

The two paths intersect in Figure 3.8 at point 0. 

The letter S represents the present position of the 
satellite and the letter T represents the present position 
of the target. These three points <0,S and T) on the 
surface of the sphere define a spherical triangle. 

Spherical trigonometry is used to attempt to find a suitable 
coordinate system for the satellite tracking analysis. 

B. COORDINATE SYSTEMS 

Several coordinate systems are tested at this time. The 
set of reference axes is a relative one. It is assumed that 
the origin of the coordinate system is moving with the 
satellite relative to the earth. Appendix B explains. 

1 • Us e Variab le s £ and 6 

Supposedly any two independent variables from 
Figure 3.9 can form the basis for a state-space model for 
target tracking. The variables 1 and 6 are chosen for the 
first attempt to form a suitable coordinate system. The 
original four state variables x(l) through x(4) are defined 
below. Derivatives are found for the four state variables. 
In an attempt to establish a least-complicated structure, 
however, these lead to other state variables and other 
derivatives ; i . e . , 
x(l) = cosine($) 
x(2) = sine($) 
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x < 3 ) 



cosine <6 ) 



x ( 4 ) = sine<6) 

Taking derivatives of two of the above variables yields 
x(l> = -x ( 2 > 
x(3) = -x<4) 

Using the trigonometric law of sines X(2) can be defined as 

x(2) = / sine (1\> sine ( t) j = sine (F\) sine(6 ) 

sme(6 ) "J sine ( <x) 

Since in this case a and sine(oc) are constants the second 
identity is easier to use. Sine(a) is replaced by the 
letter A and sine(6) is replaced by x(4) . This yields 

X (2) = sineC'rL>X<4) /A 

The same procedure is used for X(4>. The law of sines 
provides an easy alternate form. X(4) can be defined as 

x<4) = sine(T) A/sine<9) 

Table 3 summarizes the mathematics of expanding the state 
space. It is soon obvious that this state space is not 
readily reducable to a single structure by increased 
dimenssion. It therefore is abandoned. 

2 . Use of Variables 8 and S 

Two new independent variables are now used in an 
attempt to generate a state-space model. As before the sine 
and cosine functions are employed. 

X(l) = sine(5) 
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TABLE 3. 



STATE SPACE FOR VARIABLES $ and 6 



X(l> = cosine(S) X(2> = sine(iE) 

X(3> = coaine(6) X(4) = 3ine(6) 

X(l) = -X (2) and X(4> = -X<3) 

X(2) = sin(T\)sin(6 ) /sin(a) = sin (T|> X C4 ) /A 

X (2 ) = co3<r\)T\X(4)/A ♦ sin(T\)X<4) /A 
Let X < 5 ) = cos(T\) , riX<4) and X(6> = sin(T|)X(4> 

X (5 ) = -sinCrp <t\) 'XC4) + cos (T\)T\X ( 4 ) * cos <T^ <T\) (Tp X < 4 > 

Let X ( 7 ) = sin <T\) < 1\) * X ( 4 > and X<8) = cos (TV^X < 4 ) 
and X < 9 ) = cos <T\> <T\> <f[) X ( 4 ) 

X < 7 ) = cosOT!) <T\f X<4) * sin <T\) < 2I\> fjx < 4 ) + sin <T\) CT\> J X < 4 ) 

Let X<10) = cosO\) <T\> 3 X<4) and XCll) = sin <i\) < 20J>nX (4 ) 
and X ( 12) = sin (T\> CH) 1 X < 4 ) 
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X(2) 



cosine ($) 



X(3) = sine(9) 

X<4) = cosine(9) 

The derivatives for two of the above are 
X (2) = -x(l) 

X (4 ) = -x (3) 

As before sine (a) is represented by the constant A. 
The law of sines is also used again. Table 4 summarizes the 
mathematics of the expanding state space. As before these 
variables as a basis for a coordinate system do not seem 
suitable for tracking purposes. 

3 • The Proiected-Qrbit Technique 

This technique involves projecting the circle which 
is the target track onto the equator which is the 
satellite's ground track. Figure 3.10 illustrates the 
technique. The target position T is projected down onto the 

satellite's orbit at position T' . 

*• 

Let t = w(s) (the angular velocity of the satellite) 
and let $ = w(t> (the angular velocity of the target). 

For convenience let Y = ST' (see Figure 3.11) and X 
= ST' the arc length (see Figure 3.12). The arc length 0T' 

= 0T and angle 2 = angle 5' . It follows that 

5 = t ♦ 8 

and 

0 = 5 - T. 
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TABLE 4. STATE SPACE FOR VARIABLES 3 and 0 



X < 1 ) = sine < 3 ) 

X ( 3) = sine < 9 > 

X < 2 ) = -X<1> 

X(l) = sin (T)) / sin (8) 
sin <a ) 



X<2> = cosine(I) 

X(4> = cosine(8) 

X < 4 ) = -X < 3 ) 

X(3> = sin ( a) ( sin ( 3 > \ 
sin (T\_> j 



X<1) = X < 3) sin (T\) / A + cos CT|>7^X ( 3) /A 
Let X ( 5 ) = X(3)sin(T\) and X ( 6 ) = cos Cr\>T\X ( 3 ) 

X (5) = X ( 3 > s i n C1\) + cos (T^TlX < 3) 

Let X < 7 ) = X ( 3 ) sin (f\) and X(8) = cos <T\)f}X < 3 > 

X <7> = X(3)sin<n.) + X ( 3 ) cos = X(3)cos(T\)l\; <X<3> = 0) 

Let X <9> = X<3>cos<n>T\ 

• ' J • • r ■ „ 

X (9) = X(3)sincr\> (T)_) 1 + X (3) cos(T\)l\ 

Let X(10> = X(3)sin<T\) <f\) 1 and X(ll) = X < 3 ) cos (IQ) l\ 

* -■ - -i 

X(10) = X < 3 > cos <T\> CT\> + X ( 3 > s i n <T\) < 2f\) 
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Figure 3.11 The Line Segment ST’ 




Figure 3. 12 The Are ST’ 
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Taking time derivatives we find 

» * 

8 = S - t = W(t ) - (iKs) 

To obtain a state space interpretation Y and X and 
their derivatives are defined. Figure 3.13 shows an angle 
bisector which splits 0 and ST' . This forms two right 
triangles. Using geometry and trigonometry it follows that 

sine ( 0 /2 ) = T'A/R 
AT' = ST'/2 
ST' = 2Rsine (0/2) = Y 
Y = 2rcosine ( 0/2) 0/2 = Rcosine ( 0/2 ) 0 = 

Rcosine (0--r) (w(t)-w(s) ) 

X = ST' = 0R 
X = 0R > R0 f 

but R = 0; 

therefore, 

X = R0 = R ( w < t ) -w < s ) ) 

Recall R is a constant for this case and therefore R 

is zero. 

With Y and X defined above a logical state space is 
originated. The measurement equation is also important. 

In order to define a measurement equation 
Figure 3.14 is constructed. Note that points 0,S and T' do 
not form a straight line. 

0T' f OS + ST' 

However, in terms of arc lengths the following applies; 
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OT' = OS ♦ ST' 

The dotted line bisects angle a and therefore 

BT = TT'/2 

Using geometry and trigonometry it can be shown that 

TT' = 20T ' sine ( a/2 ) 

and 



BT/OT = sine(a/2) 



From Figure 3.10 it follows that 
"X OT' = 2Rsine ( <p + t ) /2) 

A new triangle is defined in Figure 3.15. Using simple 
geometry it is obvious that 

= tt/ 2 - a/2 , 

and by simple trigonometry 

DT = TT ' sine ( 6 > 

and 

DT' = TT'cosine(b) 
tangent(€) = TD/SD 

therefore 



SD = TD/tangent (€) 

The line segment ST' is solved for 

ST' = SD + DT' = TD/tangent ( € ) + TT'cosine(b) 

It is essential to express ST' in terms of usable variables. 

ST' = TT '( sine ( ft) /tangent (€ ) + cosine(D)) 

is achieved by substitution for segment length SD . 

TT' = 20T'sine(a/2) = 2(2Rsine( (ct + t)/ 2) )sine(a/2) 



therefore 
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Figure 3.13 Angle Bisector 





Fi g ure 3. 15 



New Triangle 



ST' = 4Rsine o + -r j sine(a/2) alne(fl) + coaine(fl) 

^ 2 , tangent (€) 

Let 4Rsine ( ( a + T ) /2 ) sine < ot/ 2 ) = Q for convenience. 

Then ST' = Q ^sine (15 ) /tangent ( €) * cosine(G)j. Also, 

recall that 6 = tt/ 2 - a/2 

ST' = Q | sine (ir/2-a/2) + cosine ( ic/2-a/2 ) 
tangent (€) _J 

It is important to relate X and Y. 

X = ST' = 9R and Y = ST' = 2Rsine(9/2> 

For convenience use a = X/R 

Y = 2Rsine(x/2R) and R = X/a = Y/2Rsine (a/2) 

therefore 

X = aY/2sine<a/2) 

It is now necessary to solve for € which indicates bearing. 

Y = ST' = Ocosine ( a/2 ) cotangent <€ ) + Qsine(a/2> 

cotangent<€) = Y - sj.neKa/2) 

Qcosine(a/2> cosine(a/2) 

Taking the arccotangent of each side, substituting back in 

for 0 and then simplifying leads to 

€ = arccotangent s in e (X2RJ__ - tangent (a/2) 

sine(5/2)sine(a) 

< — 

This is the measurement equation. The final result is a two- 
variable state space and a measurement equation. From the 
measurement equation an observer can be modeled. Table 5 
summarizes the results. 
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TABLE 5. STATE SPACE IN X AND 




Table 6 summarizes the results for an alternate 
coordinate system . 




A previously researched thesis [Ref 63 provides a 
useful coordinate system based on latitudes and longitudes. 
Table 3.4 summarizes this coordinate system. 



54 



IV. OBSERVABILITY 



At this time it is important to determine observability 
of the target. Assuming the satellite's detectors yield 
only the bearing information (as developed in Chapter 3), 
it is necessary to obtain the target's position defined with 
X and 5. Since this cannot be measured directly it must be 
estimated in some way. Analyzing the state equations from 
chapter three it appears that if and are given the 
knowledge of X(0) and 1(0) (the initialed states) are 
sufficient to determine X(t) and !(t) all along the 
trajectory . 



The possibility of estimating the state from the output 
measurement under the assumption that the system dynamics is 



the observability of the system must be determined. 

All changes in a given system must be reflected in the 

output for the system to be observable. Figure 4.1 shows a 

# 

simple example. The states X and X are estimated from 
knowledge of the output. If all states are observable for 
all times then the system is completely observable. This is 
defined in detail in [Ref. 6] . 




O 



completely known is defined as observability. Therefore 
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Th© test for observability on linear systems is fairly 
straight forward. In general the observability of linear 
systems is established by examining the linear 
independencies of the columns in the matrix function 
C(t)$(t,t c )# where CCt) is the measurement matrix and 1 ( t , t.) 
is the state transition matrix. For the linear time 
invariant system this leads to a simple rank test on the 
appropriate observability matrix. To illustrate this 
approach we may consider the following linear system, 

X = AX 
Y = CX 

where A is a square matrix, X is an n-dimensional state 
vector and Y is an m-dimensional measurement vector. 

Assuming that Y(t) is differentiable up to the (n-l)th 
order, we have 

Y'= CX = CAX 
Y = CAX'= CA 2 X 




A matrix is set up as shown: 
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c 





i r 


Y 




Y' 


= 


Y 1 '”' j 





CA X 



CA 



=> Y = Q X 

c 



where Y is the extended measurement vector o£ dimension mn . 

The matrix Q is the observability matrix and has 
dimension mnxn . Clearly if rank Q c = n„ X can be uniquely 

A 

solved in terms of the measurement Y. It should be noted 
that differentiability of YCt) is not required for 
observability of linear systems and that the above result 
can be derived in a different way. The presented procedure, 
however offers the possibility of extending this result to 
the observability of non-linear systems; namely in the 
equation, Y = Q 0 X, any subset of n linearly independent 
equations can be considered Y = OX , Q. can be regarded as the 
Jacobian J of Y with respect to X. Hence if rank J = n the 
linear time invarient system is observable. 



A. OBSERVABILITY OF NON-LINEAR SYSTEMS 

In the past it has been very difficult to determine 
observability of non-linear systems. Recently, however, a 
number of new methods for obsevability testing have been 
developed . 

A fairly simple test that determines observability in 
non-linear systems is developed in CRef . 6] . This new 
method tests observability and identifies any unobservable 
states that may exist. If the system is unobservable than 
the measurement doesn't provide enough information for state 
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estimation. System observability and state estimation are 
directly related. 

A non-linear system can be represented as follows: 

X(t) = fCXCt), u(t),t> (4.1) 

A r 

where f() is an n-function, XdR , UCR . 

The measurement equation is: 

Y ( t ) = h ( X ( t ) , t) (4.2) 

where h() is an m-function, Y <3 R"\ 

Y(t) must be differentiable up to the (n-1) order and 
U(t) must be differentiable up to the (n-2) order. A state 
X. (t 0 ) is observable at t if knowledge of the input U(t) 
and the output Y(t) from time to to t , enables X (t c ) to be 
determined • 

Equation 4.2 is differentiated to provide a necessary 
system of non-linear measurement equations. 



Y = h ( x ) 

y ' = 0 h + J h J x = h T «■ h x f = h, (x,u) 

^ t J> x . j t 



*< 

li 


dh' ♦ h , 


Jt 


+ -S h , 


OU 


= 






o x 


viX 


J u 








h , 1r 


+ h 1<x f 


♦ h 


= 


h. 


(x. 


U,U' > 


A '1 

Y = 


^ tiw J 


+ 


X 


t + 


^ (n 


-i| • • • 




i*-3J 




~ ^ n-i 


< X , 


XX j 


. / n-2. 

t 4 - - - - ) 



Define an mn measurement vector Y by 
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and an mn function HC) to be 



H <•> = 



h 



(\A 



The functional relation in vector form is 

Y = H ( x , v ) (4.3) 

where v(t) is a function of , i = l,....(n-2). 

With respect to equation 4.3 the question of observability 
can t?e developed as the existence of an inverse of this 
function. Clearly, this is related to the Jacobians of H. 
However, for the case of linear systems an inverse, if it 
exists, has to be unique. It will be shown late that here 
this need not be the case. 

Two conditions must be met for this non-linear system to 
be observable in the strict sense. These two conditions are 
connectedness and univalence. Connectivity is a necessary 
condition for observability in the strict sense and a 
necessary and sufficient condition for observability in the 
wide sense. Connectedness is satisfied if every state is 
connected to the output in some way. Uni valence is 
satisfied if every state is uniquely determined in terms of 
the measured output. A one-to-one mapping is considered 



univalent. Both conditions together are necessary and 
sufficient to establish observability in the strict sense. 

To explain this consider Y(t> as expanded in a Taylor 
series . 

Y ( t ) = y<t c > + y(t c ) (t-t, ) + Y*y ( t c > ( t- t c + .. 

+ ( 1/ <n-l ) ! ) y Lft < t c ) < t-t c > + r<t) 

Knowledge of the measurement trajectory Y(t) is 
equivalent to knowing the coefficients and remainder in this 
equation. This Taylor series expansion is considered unique 
and therefore the coefficients are also unique. The 
coefficients are the elements of the measurement vector Y. 
Therefore, any state (t) is observable as long as it is 
connected in a one-to-one manner to an element of Y. If the 
connection is not one-to-one the system is observable only 
in the wide sense. An example of this follows. 



X * = Y 




In this example X can equal either the positive or 
negative square root of Y. This is a multiple valued 
function. Therefore the mapping is not one-to-one. Hence 
the system cannot be observable is the strict sense. 

As mentioned earlier the connectedness condition is 
satisfied if the existence of an inverse of the function 
can be established. According to [Ref. 6] the inverse 
function is considered established if the determinant of the 
nxn Jacobian of H does not equal zero for all x and v. If 
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J = 0 then one or more states are unobservable. If H has 
More than one right inverse then univalence is not satisfied 
and the system is only observable in the wide sense. 

However this will suffice for this case. 



B. JACOBI ANS AND NON-LINEAR ANALYSIS 

Recall from chapter three that two possible bearing-only 
measurements have been found. These are illustrated in 
Tables 5 and 6. Both the measurement equations (bearing- 
only) are highly non-linear. The Jacobian matrices in 
Tables 7 and 8 are defined using the measurement equations. 

These Jacobians are necessary for observability 
analysis. It is essential to ascertain whether or not their 
determinants are equal to zero. The Jacobians in 
Tables 7 and 8 are based on the state space and measurement 
information available in Tables 5 and 6 respectively. 

The calculations for r? , X# c5^ ^ and °>v 

concerning the variable X are in Tables 9 through 11. The 
result of the calculations for , , j ^ /J y ,0^ /V>^/ and 

concerning the variable Y proved very complicated and 
are therefore not used. 

It is obvious that the determinants of these Jacobians 
are going to be difficult to analyze. Table 12 illustrates 
the determinant solution. To discover whether or not a zero 
value is achieved two separate methods are invoked. First, 
a simplification is utilized. Recall from chapter one that 
a satellite is restricted regarding the amount of area it 



61 



can cover. Therefore the numerical value of the angles $ 
and 'x are restricted. The variable X is composed of <5 and <r 
(X = R($-'t)). The development of the spherical triangle in 



TABLE 7 JACOBIAN UTILIZING X 
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TABLE S 



JACOBIAN 



X 



PART 1 



T- 
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TABLE 10 



JACOBIAN - X- PART 2 
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TABLE 11 



JACOBIAN 



X 



PART 3 




65 



TABLE 12 



DETERMINANT OF JACOBIAN 
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chapter three assumes a limit to angles 5, t and oc. For the 
purposes of this work the range of possible values for these 
angles are confined as follows; 

0 (plus epsilon) radians - tt/ 4 radians 

t: 0 (plus epsilon) radians - tc/ 4 radians 

a: 0 (plus epsilon) radians - tt/2 radians 

These limitations are reasonable for this particular 
case and they enable a useful simplification. For small 
angles the sine of that angle can be approximated as the 
angle value itself (in radians). Table 13 verifies the 
assumptions made here. 

TABLE 13 . SIMPLIFICATIONS 

tt/ 8 = 5/2 at 5's maximum 
tv/8 = .3926991 radians 

sin(ir/8) = .326834 

\ > 

For this case the angle a is held constant. Using all 
the above assumptions and approximations the Jacobians and 
their determinants are greatly reduced in complexity. 

Tables 14 and 15 summarizes the development of the 
simplified Jacobian and determinant for the measurement 
equation containing the variable X. This Jacobian is chosen 
because it is the least complicated of the two. 
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TABLE 14 



SIMPLIFIED VERSION 



Measurements : 




j 

I 


% = arccot/ X-2R5S 


S = ainMa/2) 




l R5a 


- a = sin(a) 




II 

1 

& 


/ xi - ix \ 




( 


X J - 4R5SX «■ 4R J 5 J S^ 


Jacobian : 




j 


J lt 




I 

aR ( 


\ \ 


X 1 - 


4R<5SX + 4R , 1 1 S 




n 

£ 


aR . (x i 


/o^ 


X 1 - 4R5SX ♦ 4R J $ J S ' J 


J A = O j //C = 


aR 


/ 2X1 < X-2R1S) «• 5<4R*5*S-X*j 


//J X x* 


- 4R$SX ♦ 


4R i i , SV. x 1 - 4R1SX * 4R , $*S 


= ^7/i// = 


aR 


/. > ! 

i X < 4R *5* S-X J ) +4RSX5 < X-2R5 ) 1 


x* 


- 4R5SX + 


4R 1 ! i SL X a - 4R1SX + 4R J S*S 
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TABLE 15 



LINEAR DEPENDENCE 
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The second method involves the use of computer 



simulation. The determinant (with the variable X as before) 
is calculated using no assumptions. The angles 5 and nr are 
incremented in small steps from 0 radians to tt/ 4 radians. 

The angle a (that is actually constant) has been given five 
different values for test purposes. The value of the 
determinant approaches zero but never actually achieves 
zero . 

Note from table 15 that. 

jJj = a 2 (XI - 5X) (4.4) 

X* - 4RSSX + 4R 2 $> 2 S 

If XI = <5X the determinant of the Jacobian is zero and 
the system is not observable. From Table 15 it is apoarent 
that the system is not observable if X and <£ are linearly 
dependent . 

The various calculations and computer simulations 
indicate that the system initially developed in chapter 
three (Table 5) is observable under general conditions. 
According to the results of the computer simulation the 
determinant never actually reached zero for the particular 
test values (50 test values) chosen. However there are 
certain conditions dependent upon the initial values 1 and 
't Q when the system is not observable. Equation 4.1 defines 
these circumstances. Since the system is observable for 
most cases it is possible to design an adaptive observer. 
This is the goal of Chapter 5. 
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V 



OBSERVERS 



At this point in the satellite tracking project it is 
necessary to approach the state estimation problem. In 
general, a process can be characterized by the state and 
measurement equations ; 

X(t) = AX ( t ) (5.1) 

2 ( t ) CX(t) (5.2) 

where X(t) is an n-element column vector representing the 
states, Z(t) is a q-element vector of measurements, A and C 
are nxn and qxn matrices respectively. The order of n is 
greater than or equal to the order of q. It is desired to 
obtain knowledge of the states of X(t). However, it is 2(t) 
that is measurable and therefore known. For this reason it 
is necessary to estimate X(t). The device for estimating 
the state X(t) is the observer. 

A. LINEAR LUENBERGER OBSERVERS 

A model of an open loop observer can be constructed as 
follows ; 

Y ( t ) = AY ( t ) (5.3) 

This model operates as does the process or plant. The 
estimation error is 

e ( t ) = Y ( t ) - X(t) (5.4) 
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This is differentiated with respect to time(t) . 



Therefore , 



e ( t ) 



Y ( t ) - X(t) 

AY ( t) - AX(t) 
Ae < t ) 



e(t) = e ^ e<0) 



9 (t )e(0) 



(5.5) 



(5.6) 



If the eigenvalues of A all have negative real parts 
then e(t) approaches zero as t approaches infinity. The 
error signal e(t) decays at a rate determined by the 
location of the eigenvalues of A. The rate of response of 
the process is also determined by the eigenvalues of A. 

To make the decay rate independent of the dynamic 
process, often a state model is 'driven' by an error signal. 
In this case the observer is characterized by; 

Y(t> = AY < t ) + G ( 2 < T ) - CY < t ) ) (5.7) 



and the error equation is 



e ( t ) = Y ( t ) - X(t) 



The matrix G is an nxq gain matrix which can be selected 
to determine the rate at which y(t) approaches X(t). The 
observer equation can be written as; 



Y = AY ( t ) + G ( ( X ( t ) -CY ( t ) ) = AY (t ) ♦ GCe(t) (5.8) 

The driving term is GCe(t) and its purpose is to drive 
the estimate Y towards X. The error equation is 
differentiated with respect to time. 



e ( t ) = Y ( t ) - X(t) 

e ( t ) = AY ( t ) +G ( CX ( t ) - CX(t)) - AX(t) 
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e(t> 



(A-GC) ( Y ( t ) -X(t) ) 



(5.9) 



e(t) = Fe(t) (5.10) 

where F = A-GC 

e ( t ) = e pr ' e (O ) ' (5.11) 

As expected, the results depend on the initial 
conditions. The gain matrix G must be chosen to place the 
eigenvalues of F at suitable locations. In designing an 
observer it is usual to place the eigenvalues of F where 
desired and then determine G. 

If the process has an input signal, u, then the 
appropriate equations are; 

X(t) = AX ( t ) + BU(t) (5.12) 

2 ( t ) = CX(t) 

The observer for this system is characterized by the 
following equation . 

Y ( t ) = AY(t)+BU(t)+G(Z(t)-CY(t) ) (5.13) 

The error response is the same as in the above and is 
given by equation 5.11. A block diagram of the observer is 
shown in Figure 5.1. 

When designing a linear or non-linear observer it is 
desirable that the observer error become small rapidly and 
that the observer not be very responsive to noise. It is 
very difficult to meet both these goals since they seem to 
conflict. The observer must also be supplied with a set of 
initial conditions. 
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If some of the states of a given system can be measured 



there is no reason to estimate them. An observer of reduced 
dimension can be designed to estimate only those states that 
cannot be directly measured. For example, if there are n 
states in a system and q of them can be measured, the 
observer only requires (n-q) states. Figure 5.2 illustrates 
the reduced order observer. 

B. NON-LINEAR OBSERVER 

For the design of the non-linear observer the techniques 
developed in [Ref. 7] can be applied. 

Declare two given n-dimensional vector valued functions 
of time to be Y(t) and 8<t). Let these two functions be 
related by H(0) as Y(t) = H(0). The vector H(8) is 
considered known. 

Assume that each ith row H<0) is differentiable at least 
once on all parameters 0. Therefore an nxn Jacobian matrix 
is defined as 



J ( 0 ) 




J; • (0) = J H; < 8 ) 

J . 10 - 



(5.15) 



Also ; 



Y ( t ) = Hf(8) ♦ JH(8) i8 = J(8)8(t) (5.16) 

J 0 ' J "2- 
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Figure 5.1 Basic Observer 




Figure 5.2 



Reduced Order Observer 



Figure 5.3 is the general non-linear observer referred 



to in CRef. 71. The transpose of the Jacobian is represented 



and 0^ (t) represents a solution of Y(t) = H<0) for any given 
trajectory Y^ (t) . In regards to this figure there exists a 
positive scalar s and a time T>0 such that for time t>T 
<!| 0 S (t) - 9 ^ <t)|j < s) if the following three conditions are 

met • 

1. K is positive definite 

2. The magnitude of the determinant of the Jacobian is 
bounded both from above and away from zero from 
below for the particular trajectory Yj (t) . 

3. The magnitude of Y d (t) is bounded from above. 

In addition, by increasing the minimum eigenvalue of K 
the scalar s is made arbitrarily small. 

An estimation error signal is defined as; 



as J (0). A string of parallel integrators are denoted 




e(t) = 9 3 ( t ) - 9 d ( t ) 



(5.17) 



Rearranging provides ; 



0 3 (t) 



e(t ) «■ 9 d <t) 



Referring to Figure 5.3 the following relation is 



stated ; 



0 S = -KJ 1 < 0 $ ) <H ( 0j, ) - H < 0 j > > 



(5.18) 



For notational convenience the following is defined; 



Z(e,0 d ) = H <e + 0 j ) - H(0 d ) = H<0.) - H(0^) 



(5.19) 



In this case sigma (S) does not refer to a summation. 



It is simply a variable. 




Figure 5.3 Non-Linear Observer 
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It follows from the above calculations that: 

e = 0 3 - 0^ = -KJ T (0 5 ) Ke,0 d > - 0j C5.20) 

The time-varying Lyapunov functions of the error signal 
are defined as: 

V(e,t) = KZ ‘ (e, 0^ )Z(e, 0^ ) (5.21) 

Since the determinant of the Jacobian is bounded 
Z(e,8^> = 0 if e = 0 then it follows that V(e,t> > 0 if 
e = 0. Further calculations are summarized in Table 16. 



The minimum eigenvalue of J(e + 0j)J (e+0j) is denoted by 
which is always greater than zero. The minimum 
eigenvalue of K is denoted by /\ ^ which is also greater than 
zero so both eigenvalues are positive. It follows that; 

V < -)\f /) K ]|Z<e,0 d )[| '-0, >Z<e,0 d > (5.24) 

Y d = J ( 0 j ) Qj (5.25) 

6 A = J<e J > ' Yj <5.26) 



1 -( 

Since Y and J(0 ) are both bounded J<0 d ) and 0j are 
also bounded for all t > 0. Therefore |j 0j |j ^ C 4 and 



i|J<0 d )|| < C, 

V < - .^^|Z(e,0 i )||' - C ( C,ljZ<e,0 d )jj 
V < J Z < e , 0 j >|j -C.Ci /^j 2/\f^ K - > Ci/2/lv/l^ (5.27) 

C. DEVELOPMENT OF SATELLITE OBSERVER 



For the design of the satellite observer (the system 
initially presented in Table 5> the state 0 as described 



both X 


and 


5 


0 = 


x, 


= 


X 








_ 0 _ 



(5.28) 
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TABLE 16 



LYAPUNOV FUNCTION OF THE ERROR SIGNAL 



V(t) = 



(e/t) 
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Oe 
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<0 £ 
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Cw- 

u = 6^ 

r 

Y = ( 

A 

The results developed earlier in this chapter 
(part B) are now applied to the particular problem of the 
satellite tracker. The Jacobian of Table 7 is redefined in 
equation 5.31. 

It is intuitively obvious that; 

Y = H(X,5,U> = H ( 0 ) 

The further development of both H(8> and 0 is presented 
in table 17. The information in Table 17 leads to the 
design of an observer for the satellite. This design is 
depicted in Figure 5.3. Figure 5.4 shows the model used to 
check observer results. Both the observer and model are 
simulated using the FORTRAN programming language. The 
results are discussed in the next chapter. 
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TABLE 17 



CALCULATIONS FOR OBSERVER DESIGN 
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Figure 5.4 Satellite Observer 
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VI 



SIMULATION RESULTS 



The basic observer designed in Chapter 5 is simulated to 
prove authenticity. Recall that the original state and 
measurement equations were developed in Chapter 3. Due to 
the satellite's limited GSW, several parameters (including 
angles 5 and tO are confined in their maximum attainable 
values. Therefore, these parameters were kept within 
specific ranges. The angle alpha (a) is taken to be a 
constant 30*. This value is selected for convenience. 

The angular velocity of the target (oa-> is chosen to be 
0.0005 rad/sec. This sets the target speed at about 
2 mi/sec. The angular velocity (GU;) of the satellite is 
chosen to be 0.00045 rad . sec which yields a speed of about 
1.8 mi/sec. These are arbitrary but not unrealistic values. 
A. THE BASIC CONTINUOUS OBSERVER 

Figure 6.1 depicts the basic simulation block diagram. 
The initial conditions $ and X are 0.2 radians and 1000km 
(621.5 mi) respectively. These values are arbitrary but 
within the acceptable range. Recall that X = R(3 I-t). 

Figure 6.2 is a block diagram depicting the basic 
observer design that is simulated. Note that Figure 6.2 is 
virtually identical to Figure 5.5. 
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PR0C&S5 

SATE- LLlTt 

O. r\J 

TARGET 

' l Xe, (/> 0 , CO jc. j ^Cj- <xr© 



JL 

jL 



OB5ER.V BR 



-> A 

A 

- S P 



Figure 6- 1 Basic Simulation Block Diagram 




Figure 6-2 Observer Block Diagram 
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The equations representing the various elements in 
Figure 6.2 were initially developed in Chapters 3,4 and 5. 



Table 18 summarizes these equations. 

The gain matrix (as developed in Table 17) is chosen to 
be a diagonal matrix of the form: 

K = (k, 0 

0 

The elements k and k both must be greater than zero. 

A diagonal matrix was chosen because it is fairly easy to 
work with and it allows the error changes to be monitored. 

There are two ways to describe the observer to be 
simulated. It can be described in the continuous domain or 
in the discrete domain. 

The observer was initially represented as a continuous 
model. A forth-order Runge-Kutta method is used for 
integration. The continuous model was used in order to 
eliminate any possible errors that could be introduced as a 
result of discretization. However, simulations for both 
versions produce no real differences. 

B. THE DISCRETE OBSERVER 



With the advent of micro-computers a discrete observer 
acting as a one-step predictor is a more realistic choice. 
In this model the foreward difference approximation is 
employed . 

The basic definition is; 

X(t , ) = X(t K ) +4t X(t c ) 

l<-r l ^ s 
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TABLE 18 



OBSERVER EQUATIONS 
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1 



Full Non-Linear Observer 



First the original, highly non-linear observer 
(figure 5.5) was simulated using a FORTRAN program on the 
IBM 370 system. All attempts to find values of gains and 
integration step for which observed states converge were 
unsuccessful. Typically the observer enters steady state 
with very small errors e, and e^; however, the estimated 
states were far from actual values. It is possible that 
this behavior is the consequence of the sinusoidal periodic 
functions (sin(x/2R) and sin($/2)) in 7 ? and ^ . Recall from 
Chapter 4 that this observer is not shown to be observable 
in the strict sense. Univalence may not exist due to the 
reiterative nature of the sine function. In order to avoid 
periodicity an approximation of a non-linear observer is 
considered . 

The simplified observer has the same form as the 
full observer (see Table 18) , but the measurement equations 
and the corresponding Jacobian are approximated as developed 
in Chapter 4 (and repeated in Table 19) • 

The simplified observer simulated was being very 
insensitive to changes in X and 5. The observer stabilized 

A A- 

with a small error (e ( and e a) while X and <5 remained 
significantly different from X and $ respectively. The 
exact reason for this behavior is unknown. However it is 
theorized that the relative smoothness of the arccotangent 
curve (arccot^?) is responsible for this behavior. 
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TABLE 19 



SIMPLIFIED OBSERVER 
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In order to circumvent this problem the measurement 
signal is redefined. 

2 • Redefined Non-Linear Model 

A new measurement variable is defined as: 

/^ u (t) = x(t)/$(t) = Rsin(a) (cot(^(t) ) + tan(a/2) ) 

If the value of a is known and ^ (t) is measured; 
then u(t) can be calculated. Taking into account the 
assumed range for $(t) and X(t) and the values for 
parameters R and a, the expression X/ (R$sin(a) ) -tan (a/2) 
never achieves a very large value. Therefore <t) is 
bounded away from zero. This implies that cotan(^(t) ) has a 
bounded value. 

The derivative of^-iXt) is: 

^ (t) = Rsin (a) <-^(t)/sin J <^(t) ) ) 

where ^(t) and ^ ( t ) are measured values and y<^(t) is 
evaluated . 

The redefined non-linear model and observer are 
presented in Figure 6.3. It is noted that this observer has 
the same form as the observer given on the block diagram in 
Figure 6.2. Observer equations are shown in Table 20. 

All attempts to find the power K t , K 4 and At (the 
time increment) that would force the states to converge 
failed. A satisfactory theoretical explanation for this 
behavior was not established. However, one possible 
explanation is that e and e are extremely sensitive to 
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TABLE 20 



MEASUREMENT EQUATIONS - REDEFINED NON-LINEAR MODEL 
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change in states which results in an extremely slow 
convergence with a highly oscillatory transient response. 

3 • Redefinition of State Variables 

Analysis of the model with redefined measurements 
reveals that the form of the equations allow redefinition of 
state variables. This may reduce the sensitivity of the 
model . 

A new state vector CU^VD^is defined as: 

U = X/5 and V = 1/5 

Table 21 illustrates further development. 

The block diagram in figure 6.4 shows the discrete 
observer. A simulation of the discrete observer was run 
with = 0.0005 and =0 . 00045 . Several different gains 

were used to test convergence. Convergence is obtained for 
the following ranges: 

0 < K, At < 2 and 0 < K^At/|J z aj < 2 

A /*• 

Arbitrary initial values for U and V were employed. 

A * a. 

It seemed natural to assume Uo = 0, X =0 and V = 1 <5 =1 

is within the allowable value for $> . 

To examine the observer's ability to track varying 
angular speeds it is assumed that and are sine waves. 

= A sin(Dt) 

OJs = 3 sin(Dt> 

Different amplitudes and periods are tested and 
successful tracking is achieved. Appendix C contains the 
resulting computer graphs. 
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1 




F igure 6.3 Redefined Model 



n- 





Figure 6. A Discrete Observer 



92 



TABLE 21 



OBSERVABILITY CONDITIONS AND DISCRETIZED OBSERVER 



Measurement equations: 

/x<U, V) = h | CU, V) = U 
AJU U,V> = h^OJ,V> = V(X - $U> 

Jacobian : 



Cl 

II 

II 




J,, = ^ h, = 


0 


ju 




d V 




Jz. = i h- = -$V ; 




J i2 = dh< = X 


- $u 


JU 




JV 




Observability conditions 


are the same: 




jjj = X -$U * 0 


X 


- IX f 0 J(, 


/ j. 


1 i 




^ X 


0 


Continuous observer : 








U = K, J„ 


e 


, + K,J, ( e. 




V = K.J U 


e 


j + K 2 J 1 Z e 2 




The discretized observer 


becomes ; 




e,K = yU (U < ' V k> - 


h 


1 ,Vp> =yu^- 


A 

Uc 


) - h.CU 


* » 


v\> = far v r <x 


- 5 



Up,., 3 U < K, dte lfc K,at • Jii<UK,V P )e zl 






VKf| = + Ki & t J zi ( U* 0 V j^) e 



2. 



S3 



c> 



As expected the amplitude of the sine wave 
representing state <A/D;B/D) which can be tracked is 
inversely related to the period of oscillation. 

4 • Model Decoupling 

Analysis of the error equations provide; 
e, ■ />< U,V> - h,< U , V ) = U-U 
e , = yUOJ, V) - hi(U, V) = V < X-5U > - V < X-5U ) 

= X CV-V) -$<UV-UV> 

The block diagram in Figure 6.4 indicates the 
possibility of decoupling the observer into two first order 
models by letting equal zero (J,* is already zero) . 

Figure 6.5 shows a block diagram of the discrete 
subsystem. The discrete system of Figure 6.6 is certainly 
stable if it's eigenvalue lies within the unit circle. 

j 1 -a | < 1 0 < a < 2 >> 0 < K, /it < 2 

Under this condition the system has a steady-state 
error of zero for a constant input u. 

Hence the system acts as positioning servo-mechanism 
and is capable of tracking the relatively slow varying 
inputs with very small errors. 

5 . Subsystem Two 

Assuming that subsystem 1 is tracking u ( u(t) = 
u<t)) the error e can be approximated as: 

(X - 5U> CV - V) 

The observer is reduced to the gradient algorithm 
where the gradient is: 
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SuB^yjTtH 1 = K, e, = K,u-K,d 

| AJ. 



CORA tSPOM Di/VG DtSCRETt riODEL 
-Ct K+1 =(|-K,At) i U K +K 1 4tU K 

= Ci-cl)cXk f aa K = 

lO- ° K, At) 




Figure 6.5 Subsystem 1 




Figure 6.6 Discretized Block Diagram 
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) = X - $U = J 

Figure 6.7 illustrates the observer. This result is 
also in accordance with the method proposed in [Ref. 63. 

Defining <X-$U)V as a new measurement the 

corresponding function h is h=(X-£U)V. The Jacobian is 

j = 6u/Jv= x-iu. 

In order to simplify the stability analysis the 
Jacobian, J, is replaced with the signum function, 
signum(J) . This is sufficient to guarantee the negative 
feedback of subsystem 2. Subsystem 2 is analogous to 
subsystem 1 and will be stable if 0 < K L At < 2. 
Consequently this system also performs as positioning servo- 
mechanism and will track relatively slow varying inputs V(t) 
with very small error. 

The presented analysis leads to a definition of the 
decoupled observer. This observer is depicted in Table 22. 
Appendix D contains the computer simulation results. 

Figure 6.8 is essentially the same as Figure 6.4 
with J^i = 0 and J x ^ = signum(J 2 i ). Therefore it is 

possible to consider this algorithm a simplification of the 
algorithm developed in [Ref. 73. However, it should be 
pointed out that the decoupled observer can be derived 
directly through analysis of the two linear subsystems 
presented above. 
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TABLE 22 



DISCRETIZATION 



= . M - U 



/N. A. 

h , <U, V) 



- V(X - $U> = 



* - h 2 (U,V) 



U = K,e f 



V = K ,Sygn(Ji 2 >e , 



After discretization this becomes 



U K., = u k * K .4te lK 



= V H * K 2 At sygn<J i 2 )e iK 



I 

I 
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Figure 6.7 Subsystem S 







Figure &. Q Decoupled Observer 
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In the decoupled model it is easier to control gains 
since the J element does not affect stability. It is 
noted that slightly more oscillatory behavior is present 
at the beginning of the simulation. This is not unusual for 
gradient type algorithms. 

Satisfactory results for the product K&t = 0.1 are 
expected if the observer is analyzed from the sampling 
theorem point of view. This presupposes that both 
continuous subsystems posses only one time constant ^ = 1/K. 
The sampling theorem requires that t < T ' / 2 and AtK < 0.5. 

Experience has shown that the best results for the 
one-step integration are obtained if the integration step 
(sampling time) is ten times less that the smaller time 
constant in the system. This means At * 0.1 (1/K) and 
AtK * 0.1. 

The sampling theorem is violated if AtK=l. As a 
consequence the observer requires much more time to extract 
sufficient information from the measurement and to start 
state tracking. The simulation results demonstrate that 
when the product AtK equals unity convergence starts after 
fifty seconds. This is approximately ten times slower from 
previous experiments where the sampling theorem is 
satisf ied . 

In the case of itK = 0.01 the system pole is very 
close to the unit circle. Therefore, convergence is 
relatively slow. 
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As demonstrated by the various computer generated 
graphs, the simulation results for this system are 
satisfactory. Therefore the selected coordinate system is 
proven acceptable . 



100 



VII. CONCLUSIONS AND RECOMMENDATIONS 



The objective of this thesis is twofold. The first part 
of the thesis involves finding a suitable coordinate system 
in which to establish a satellite tracking model. The 
second part of the thesis uses this coordinate system to 
obtain an observability analysis and design an observer. 

This latter portion is the main thrust of the thesis. 

The satellite model is based on the coordinate system 
developed in chapter three. The model is simplified yet not 
unrealistic. The main limiting assumption confines Doth the 
satellite and target to a constant speed and heading. An 
unchanging direction is important due to the desire to 
maintain a constant angle a. These restrictions imposed 

on the satellite are explained in detail in chapter three. 
Satellites and targets (aircraft) would normally travel a 
great circle path which requires a constant bearing. 
Therefore, the assumption is valid. This basic premise 
resulted in the spherical triangle that forms the basis of 
the system presented in Table 5. Once the basic model 
became established an observability analysis was 
accomplished in chapter four. The basic observability 
theories developed in CRef. 6] (for mechanical springs) has 
been successfully applied to satellite observability. The 
system proved to be observable in the wide sense. 
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Further studies in the area of satellite tracking and 
observability are highly recommended. A more realistic 
satellite model could be devised in these studies. To 
improve the accuracy of the satellite motion equations the 
earth must no longer be assumed spherical and the 
satellite's orbit must no longer be confined to simple 
circular orbit. This leads to a different coordinate system 
by necessity. Relative coordinate systems should be 
explored in more detail. Appendix B discusses this briefly. 
A possible approach involves using the Euler equations of 
motion. A target position with respect to a satellite prime 
axis -could be transformed via pitch, yaw and roll to a 
satellite normal axis. Enhanced satellite tracking studies 
should include a target that is not confined to a constant 
heading. The target should be allowed to alter course to 
further complicate the coordinate system and observer model. 
Both the satellite and target should be subject to random 
disturbances. This leads to the development of stochastic 
models. References 6 through 8 are highly recommended to 
anyone pursuing this course of research. 

This thesis has successfully developed a model for 
satellite tracking. The system has been determined to be 
observable and an observer has been designed. The observer 
model has been successfully simulated. 
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APPENDIX A 



ELEMENTS OF ORBIT 



Satellites travel in elliptical orbits. Figure A.l 
depicts a basic ellipse. The distances a,b and c are 
related by equation A.l. 

c 2 = a 2 - b 2 CA. 1) 



In this equation: a = semi-major axis 

b = semi-minor axis 
c = distance between foci 
Eccentricity is also a parameter often used in 
connection with elliptical orbits. Eccentricity, e, is a 
meassure of the elongation of the ellipse. Equations A. 2 
and A. 3 relate eccentricity to previously defined 
parameters . 



e 2 = 1 - (b/a) C A . 2 ) 

c = ae ( A . 3 ) 

Other elements include apogee and perogee . Apogee is 
that point on the orbital ellipse where the satellite is 
furthest from the earth (focal point). The perigee is that 
point where the satellite passes closest to the earth. 

These new parameters are related to a and e by equations A. 4 
and A . 5 . 

apogee distance = d * = a(l*e) (A. 4) 

perogee distance = d^ = a(l-e) (A. 5) 
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Any two of the above six parameters can define an 
ellipse. For the special case of a circle; e = 0, d = d 
c = 0 and a = b . 

To describe satellite motion more information is needed. 
The inclination angle,, i d , is the angle of intersection 
between the orbital plane and the earth's equatorial plane. 
Figure A. 2 illustrates. Information on the period of 
revolution of the satellite is needed as well as knowledge 
of the precise time that the satellite passes it's apogee or 
perogee . 

Some references use information on the ascending node to 
determine satellite paths and positions. Figures A. 3 and 
A. 4 illustrate the relationships. [Ref. 5j uses the 
parameters in table A.l as the necessary elements of orbit. 

TABLE A.l 

T - period (in minutes) w - argument of perogee 

i - angle of inclination e - eccentricity 

)\ SSP longitude at perogee t - time at perogee 



The elements listed in table A.l completely describe the 
theory of satellite motion. Perturbations are ignored. In 
practice the motion of a satellite is very nearly described 
by these elements of orbit. 
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Figure A. 2 Orbital Plane 
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Figure A. 3 Elliptical Orbit Diagram 
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Figure A. 4 Ground Tracks and Ascending Node 
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APPENDIX B 



RELATIVE MOTION 

The general theory of relative motion is developed using 
simple frames of reference. An object in frame A moves with 
respect to that frame. Frame A moves with respect to the 
fixed frame B. These frames are considered to be coordinate 
systems or sets of reference axes. Figure B.l depicts an 
object in the prime system. The X-Y axes represent the 
fixed system. In figure B.2 the prime system moves with 
respect to the fixed system a distance d. The object in the 
prime system moves a distance d' within its system. 

Therefore the object moves a distance d rcT . A l. = d + d' with a 
velocity V TcraL = V + V' . 

Further studies concerned with satellite tracking should 
involve relative motion. The fixed coordinate sytem is a 
three dimensional system with its origin at the center of 
the earth. The center of mass of the satellite is the 
origin for the relative or prime coordinate system. The 
prime system orbits the fixed system. 

The target is detected by the satellite and its motion 
is therefore meassured with respect to the satellite's prime 
coordinate system. These relative coordinates can be 
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related to the fixed coordinate system by a suitable matri 



transformation . 
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Figure B. 1 Two Frames of Reference 
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APPENDIX C 



ORIGINAL MODEL SIMULATION RESULTS 

This appendix contains the simulation for the original 
model. These graphs depict the states X and Figure 6.4 

presents the block diagram of the model simulated. The 
simulations are run for 10 seconds and 500 seconds. This 
insures that a variety of suitable graphs are available for 
analysis. In order to accompany test observer behavior 
several different gains and time increments are used in the 
simulations. Specific combinations of gain and time 
increment produce excellent results (as explained in 
chapter 6) . 

The observer model's ability to track a constant is 
recorded on eight graphs. The first two graphs (pages 112- 
113) depict state X for constant angular velocities and 
different gain values. The next two graphs (pages 114-115) 
depict state <£ for the same gains mentioned above. A gain 
(k) of ten is clearly beneficial for both X and since 
then /itk = 0.1. 

The following two graphs (pages 116-117) show two 
different time increments for state X, constant gain and 
constant angular velocities. The graphs on pages 118 and 
119 display the same two time increments for state <S. As 
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expected the smaller of the two time increments (0.01 
seconds) gives a slightly better performance. 

The last six graphs exhibit the model's time response to 
sinusoidal angular velocities. These are the results of the 
test used to determine the model's ability to track time 
varying functions. The first three of these graphs (pages 
120-122) illustrate a target angular velocity of 
.OOl-cos(Dt) and a satellite angular velocity of 
. 0009 • cos ( Dt ) . The variable D is varied in order to test 
observer response to different rates of change in angular 
velocities. The final graphs (pages 123-125) depict X and £ 
with 'a target angular velocity of .005*cos(Dt) and a 
satellite angular velocity of . 00045 • cos ( Dt) . The variable 
D is 0.01 in this simulation. This results in sinusoids of 
smaller amplitudes and longer periods. The results of 
tracking the sine wave are very sat isf actor y . It is noted 
that the relative error for both states is very small. 
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APPENDIX D 



DECOUPLED MODEL SIMULATION RESULTS 

These twenty graphs resulted from the decoupled model 
simulations Cas described in chapter S). Figure 6.8 
presents a block diagram of the decoupled model. The first 
five graphs (pages 128-132) illustrate the time response for 
state X. The time increment, a t, is held constant at 0.01 
seconds. Several different gains were used on each graph. 
The state X was followed closely by X. Different scales and 
simulation times are plotted to produce a variety of 
suitable graphs. A gain of ten consistently produces very 
good results in all five graphs. 

The next group of five computer graphs (pages 133-137) 
exhibit state 5. As above, the time increment. 7t , is kept 
at 0.01 while the gain (k) is varied. The time response for 
the lower gains in particular are very satisfactory. The 
state £ is closely approximated by <§ . 

The third set of five graphs (pages 138-142) represent 
state X. This group of computer graphs holds the gain 
constant at one and allows the time increment to vary. As 
expected, smaller time increments (0.1 and 0.01) result in 
superior tracking capabilities. 

The last set of five graphs (pages 143-147) show <£ and 
$. The gain is kept at a constant value equal to one and 
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the time increment is varied- As before, the smaller time 
increments produce the best results. 

In all of the above examples, the approximated state (X 
or $) closely tracks the actual state (X or 1). These 
simulation results are quite satisfactory. 
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